Off-shell effects on particle production* 

M. Effenberger and U. Mosel 

Institut fiir Theoretische Physik, Universitiit Giessen 

D-35392 Giessen, Germany 

UGI-99-11 

'd\ 

2i '■ Abstract 

We investigate the observable effects of off-shell propagation of nucleons in 
OO ' heavy-ion collisions at SIS energies. Within a semi-classical BUU transport 

^^ ' model we find a strong enhancement of subthreshold particle production when 

off-shell nucleons are propagated. 
> 

in 

OO 



During the past 15 years semi-classical transport models have extensively and quite success- 
O ■ fully been applied to the description of heavy-ion collisions from SIS to SPS energies (see 

O 



e.g. Refs. |ll|,0]). AH these models are based on an on-shell approximation for the nucleons 
^ I whereas the nucleon resonances are propagated with their spectral functions. However, al- 

0^ \ ready at SIS energies the collision rates are so high that this approximation for the nucleons 

is hard to justify. 

While the general formalism to treat also the off-shell transport of particles has been 

^ , available for quite some while [0-|^ only recently the interest in the in-medium properties 

pH I of the p-meson has triggered some new efforts to develop a transport theory that includes 

dynamical spectral functions (see e.g. 0-^]). In Ref. we have already presented a model 

'^ ' in which the spectral functions of the p and uj mesons are dynamically incorporated and 

'^ ■ have discussed in detail the transport equation that we solve. It is the purpose of this letter 

^ to discuss the effects of coUisional broadening also on the spectral functions of the nucleons 

and thus to treat them on the same footing as the resonances. We also want to report first 

results on particle production in heavy-ion collisions that were obtained by an inclusion of 

the nucleonic spectral function. 

Since in our method a 'spectral phase space distribution' plays a central role we specify its 
precise meaning by starting from the non-relativistic, semi-classical Kadanoff-Baym equation 



[po - H^f, g<] + [Reg, S<] = -S>(7< + T.<g> , (1) 

where the Poisson bracket is given as: 
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[x,yi = ^^-^^ , (2) 

with X and p being the spatial and momentum coordinates. The mean field Hamilton 
function Hmf for the non-relativistic case reads: 

where Mq denotes the mass of the particle and S(x,p) its retarded self energy that is related 
to the self energies S^ and S"^ which appear in the 'collision term' on the rhs of Eq. (|l]) via: 

ImS = i(S>TS<) (4) 

with the upper sign for bosons and the lower sign for fermions, respectively. 

g^{x,p) and g-^{x,p) are the Wigner transforms of the real time Greens functions and 
g^ can be interpreted as a generalized phase space density, i.e. the density of particles 
with 4-momentum p at space time point x. The Wigner transform of the retarded Greens 
function g{x,p) is given by: 

1 1 

Img = -{g> T g"") = ^a , (5) 

where a{x,p) is the spectral function, g'^ and g^ can be expressed through the spectral 
function a and a phase space distribution function f{x,p): 

g" = af (6) 

g> = a{l±f). (7) 

For stable particles / reduces to the usual phase space density. From Eq.(|T]) and the corre- 
sponding equation for g^ one obtains an algebraic solution for 5' 0: 

(8) 



The spectral function is then given as: 



bo-og--ReS)2 + i(S>TS<)2' 



(9) 



For a thermally equilibrated system / reduces to the usual Fermi-Dirac or Bose-Einstein 
distribution function; the non-trivial spectral information is then contained in a. 

In the following we will neglect the second Poisson bracket of Re^' and S*^ on the Ihs of 
Eq. (|l]). Its infiuence on the description of heavy- ion collisions is hard to estimate because 
a numerical implementation of its effect is not available so far. 

Instead of the non-relativistic Hamilton function Eq. (^ we use, in our specific numerical 



realization 110], the following one: 



H„^f = ^J{^l + sy + f , (10) 



where S{x,p) is an effective scalar potential and /x is a mass parameter that is related to 
the energy po via: 

Po = Hmf- (11) 

In the following we will express all quantities as function of /u instead of po- We, therefore, 
instead of g'^ use a spectral distribution function F{f, t,p, fi) defined by: 

Fdl.^9<'^ (12) 

and use instead of the spectral function from Eq. (|^) the following relativistic one: 
./-, - ^ _ 2 /i^rtof(f,t,p,^) 

71 [^i^ - M^Y + ^i^vi^tij, t,p, /i) 

where the total width in the rest frame of the particle Ttot is given as: 

r,„, = ^(S>TS<). (14) 

1^ 



The transport equation then reads: 

)F = S<(l±/)^-S>F. (15) 



d dHmf d dHmf d 



dt dp df df dp 
The collision term is given as: 

1 [ d^p2d^2 d^Psd/J-s d^p^d^^ 



S<(l±/)^-S>F = — |. 



2E J (27r)32^2 (27r)32^3 {2ttY2E^ 
x(27r)V^)(p + P2 -P3 -P4) lA^l' {FsF.il ± 72)^2(1 ± f)A- 

-FF2(1±/3)^3(1±/4)A} , (16) 

where Fj = F{f,t,pi, fii) and Ai, fi are denoted analogously. 

We solve the transport equation (0) by a so-called test particle method, i.e. the distri- 
bution function F is discretized in the following way: 

F(f,t,p,/i)oc5]<5(f-r;(t))5(p-p;(t))(5(^-^,(t)) , (17) 

i 

where fi(t), Pi{t), and /Xj(t) denote the spatial coordinate, the momentum and the mass of 
the test particle i, respectively. By inserting this ansatz into the Ihs of Eq. (|T3p one obtains 
the usual equations of motion for the test particles between collisions: 

(18) 

(19) 

(20) 
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dpi 
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dt 
dni 

dt 


df 
= 0. 



A formal generalization of the transport equation to a system with different particle 
species is straightforward. One obtains a transport equation for each particle species that 
is coupled to all others via the mean field potential and the collision term. In our transport 
model all baryonic resonances up to about 2 GeV as well as all low lying meson states are 
explicitly propagated. For a detailed description of the model we refer to Refs. P,p!0[]. The 
actual collision term for the nucleons is thus much more general than the one in Eq. (|16D and 
a numerical solution is only possible by making some approximations that we will describe 
in the following. 

We neglect any medium modifications of total cross sections and decay widths except of 
Pauli blocking of the outgoing nucleons. The cross sections are calculated as function of the 
'free' invariant energy y/st.^^^ which is given as: 



VSfree = V ^1 + Pcm + V ^2 + PL > (21) 

where /zi and /i2 are the masses of the incoming particles and Pcm denotes their cms 3- 
momentum. For invariant energies below 2Mjv we assume the cross section to be constant. 
We have checked that our results do not change when we use the elastic nucleon-nucleon 
cross section as function of y/sr^^^ — /ii — /i2 instead of as function of y/Sfree- The in-medium 
spectral functions of the produced particles enter only in the determination of the 4-momenta 
of the outgoing particles. For 2 particles in the final state the masses and the production 
angle are chosen according to 

;^ oc>y(9)M.-4. . (22) 

where W{9) denotes the respective angular distribution and pf is the cms final state mo- 
mentum of the outgoing particles with masses ^3 and ^4. 

Because in an actual simulation it is practically impossible to obtain statistics for a 
determination of the collision rate in 8-dimensional phase space we calculate the collision 
width of the nucleons that enters the spectral functions in Eq. (|22| ) only as function of total 
momentum in the local rest frame pirf, density p and temperature T: 

^Nir,t,p,fl) ^ TN{pirf,p{f{t)),T{f{t))). 

The temperature is determined by assuming a Fermi-Dirac distribution from the local < 
j9^ > of the test particles. The width is calculated with the nucleon-nucleon cross sections 
described in Ref. [|T0[ where we take into account Pauli blocking for the elastic part (see e.g. 
Ref . [^ for details on the calculation of the collisional widths) . In Fig. |l] we show the width 
of the nucleons as function of momentum and temperature at density 2pQ. One sees that the 
width is of the order of 300 MeV for a temperature of 100 MeV and a momentum of 1 GeV 
and therefore non-negligible. We have checked that the actual collision rates (time integrated 
and averaged over momentum) in the simulation agree very well with these widths. This 
holds also for the early, non-equilibrium stages of the collision because the collision rates 
depend mainly on the density and the momentum. We also find that it is essential to keep 
the rather strong momentum dependence of the widths, depicted in Fig. |I|. 

In Ref. we have discussed in detail that the transport equation Eq. (p!5D does not give 
the correct asymptotic solutions for particles that are stable in vacuum because a collision 



broadened particle does not automatically lose its collisional width when being propagated 
out of the nuclear environment. In Ref. [§] we have introduced a simple scheme to restore 
the asymptotic width by adding a scalar potential to the equations of motion for the test 
particles that shifts a particle on its mass shell when it propagates to the vacuum. In order 
to fulfill energy conservation in case of a system of off-shell nucleons we have to take into 
account the back coupling of the potential on the other nucleons that create the potential 
via the density. For a testparticle i this restoring potential is defined as: 

where tcr denotes the time when the particle is produced; for nucleons we have n^ac = M^. 
The mass //j is then given as: 

fJ'i (i) = fJ-vac + Si{Xi(t),t) + s{xi{t),t). (24) 

In addition to the scheme in Ref. ^ we have included the back coupling potential s which 
is given as: 

-. ^ ,x _ 1 / d^pd^ifi - fivac)F{f, t, p, /i) 

'^'' ^ 2 Jd^pd^,F{f,t,p,^,) • ^''^ 

This prescription can be formulated in terms of a modification of the equations of motion 
for the testparticles: 

^/^i h^i H'vac "^ "P ^^ H'i H'vac ^ '^J- coll fO(i\ 

dt p dt Tcoii dt 

where in the last step we have assumed that the collisional width Tcou(t) oc p{t). Since Sj 
and s enter the equations of motion as usual potentials we also get a respective modification 
ofEq. dii: 

dpi p + S fdS ^ p- p^ac - s dTcoii ds\ , . 



dt vo \ dr Trnii dr df 



This procedure conserves total energy since s is the local average of the Sj's. Numerically we 
obtain energy conservation on a level better than 1% whereas a neglect of the back coupling 
potential s leads to a violation of total energy conservation by about 10%. However, for the 
calculations of particle yields presented here the inclusion of the restoring potential Eq. (^3[) 
is only of minor importance as will be shown below. Because of this it is not essential to 
what extent our prescription mimicks the effect of the Poisson bracket that we neglected 
in Eq. (|1]). We note that the transport equation (|I]) is based on a first order gradient 
expansion and therefore does not necessarily need to give the correct asymptotic solutions if 
the gradients are too large. If the gradients are low enough already the collision term shifts 
the particles on mass shell as we discussed in Ref. p. 

In the following we exemplarily present results of the present formalism for central 
Au+Au collisions at 1 AGeV bombarding energy. We have checked that the observed effects 
are similar for larger impact parameters and smaller systems. The results that we show here 



were obtained with an EOS that has a compressibihty K of 380 MeV. With a softer EOS 
(^^=215 MeV) the effects of off-shell nucleons are quantitatively the same. In these calcula- 
tions we employ a (somewhat arbitrary) lower mass cut-off of 400 MeV to avoid unphysical 
test particle velocities. The results are stable when this cutoff is varied. 

In Fig. 1^ we show the mass differential spectrum of the nucleons for different times of 
the reaction: 

dN ^^^ f d?p 3 



where we only count the nucleons that suffered at least from one collision since the others 
contribute to a 5-function at the pole mass. The upper part of the figure displays the result 
without the potential Eq. (pS]). One sees that the nucleon distribution at the end of the 
collision {t = 40 fm) is quite strongly peaked around the nucleon pole mass (FWHM?«10 
MeV). However, there is still some quite broad tail of low mass nucleons. The inclusion 
of the potential Eq. (|23|) 'cures' this problem (lower part of Fig. ^. One should note here 
that the spectral distribution in the high density phase of the reaction is hardly affected by 
the potential. The distribution has a width in the order of 200 MeV which again stresses 
that an on-shell approximation for the nucleons is hard to justify. It is an essential property 
of these spectral distributions that they are strongly asymmetric, with an increase towards 
lower masses. This stems mainly from available phase space which gets larger for smaller 
masses even without the factor pf in Eq. (^2]) . We stress again that this is a consequence of 
the collisions and not of the potential (p3|). 

In Fig. H (upper part) we show that the off-shell propagation has a quite large influence 
on the pion spectrum. The number of high energy pions is increased by about a factor 2. 
There are essentially two reasons for this result. First, in the initial stages of the collision 
high mass nucleons can be produced which contribute to high energy secondary collisions. 



This is the effect which has been considered in Ref. |TT|. Secondly, outgoing nucleons can 
have low masses so that effectively the 'temperature' increases. This allows to produce 
heavier baryonic resonances in A^A^ —>■ NR and more energetic pions in R —>■ Nn. We note 



that also in the quantummechanical off-shell investigation of Ref. |T2[ a higher particle yield 
is obtained when in- medium correlations are taken into account. 

In the calculation of these results we have used the free A spectral function; in the 
population the correct phase space factor has been taken into account. We have ascertained 



that using the in- medium widths of Ref. |13] in the population of the A made no significant 
difference. 

The lower part of Fig. ^ shows our result for the X^-spectrum where we see that the 
off-shell propagation gives an enhancement by about a factor 2. This enhancement is solely 
caused by the increase of high energy baryon-baryon and pion-baryon collisions since we do 
not modify the total production cross sections. This would even lead to a further increase 
of kaon production because in BB —>■ NYK a finite width of the outgoing nucleon lowers 
effectively the production threshold. 

The potential Eq. (|23|) (with s=0) has only a small effect both on the pion and on the kaon 
spectrum as can be seen by comparing the dashed and the dotted curves in Fig. 0. This is so 
because in the high density phase, when the high energy mesons are produced, the collision 
rates for the nucleons are so large that the nucleons do not travel through a relevant density 

6 



gradient for which the potential would become important. The potential mainly shifts the 
low mass tail of nucleons on mass shell during later stages of the collision. For subthreshold 
particle production the exact mechanism how particles in heavy collision come to their mass 
shell plays only a minor role. The calculation which includes the backcoupling potential 
Eq. (^) and thus conserves total energy gives practically the same result (dash-dotted 
curves in Fig. ||). 

We shortly want to mention a few possible improvements for future work. By using vac- 
uum cross sections and vacuum decay widths, as it is always done in transport calculations, 
detailed balance is violated when off-shell nucleons are propagated. This could only be cured 
by working with amplitudes instead of parameterizations of cross sections and decay widths 
which would enhance the numerical effort dramatically. However, we do not expect the 
violation of detailed balance to be significant for our present results because the harder par- 
ticles for which the increase due to off-shell effects is most pronounced (Fig. |^) are produced 
in the first hard collisions. Furthermore, in all transport calculations the analyticity of the 
self energy S is neglected and the mean field potential that enters into the drift term is 
determined independently of the self energies S*- and S^ in the collision term. For off-shell 
transport this analyticity might be important for the mechanism how the particles come 
to their mass shell when they travel to the vacuum. In principle, we should also take into 
account a dynamical spectral function for the pion since its collision rates are even larger 
than the ones of the nucleons. 

In summary, we have presented a first calculation of nucleonic spectral functions in a 
transport theoretical framework and have exploited their effects on particle production in 
heavy-ion collisions at SIS energies. We have found large effects compared to a calculation 
with an on-shell approximation for the nucleons. We thus hope that this work will stimulate 
future, more refined work on the influence of off-shell properties on particle production yields 
in heavy-ion collisions. 

The authors would like to thank C Greiner, J. Knoll, V. Koch, A. Larionov, and S. Leu- 
pold for discussions and helpful comments on the manuscript. They also acknowledge a 
critical reading of the manuscript and lively arguments with W. Gassing. 
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FIG. 1. The nucleoli width as function of momentum p and temperature T at density p = 2po- 
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FIG. 2. Spectral distribution of nucleons for Au+Au at 1 AGeV. In the fower part the restoring 



potential from Eq. (23) was included (without the backcoupling potential s). 
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FIG. 3. Effect of off-shell propagation of nucleons on pion (upper part) and kaon (lower part) 
spectra. For the dotted lines the restoring potential from Eq. ( p3|) was included. For the dash-dotted 
lines the energy-conserving backcoupling potential Eq. (pq) was additionally taken into account. 
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